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The low energy behaviour of the two-dimensional antiferromagnetic Heisenberg model is studied 
in the sector with total spins S = 0, 1,2 by means of a renormalization group procedure, which 
generates a recursion formula for the interaction matrix Ag™ +1 ' of 4 neighbouring "n clusters" of 
size 2" x 2", n — 1, 2, 3, . . . from the corresponding quantities Ag . Conservation of total spin S 
is implemented explicitly and plays an important role. It is shown, how the ground state energies 



E 



(n + l) 



S — 0,1,2 approach each other for increasing n, i.e. system size. The most relevant 



couplings in the interaction matrices are generated by the transitions {S',m';n + l\Sg \S, m; n + 1} 
between the ground states \S, m; n + 1) (m = —S, . . . , S) on an (n + l)-cluster of size 2 n+1 x 2" 
mediated by the staggered spin operator Sq . 



PACS numbers: 71.10.Fd,71.27.+a,75.10.-b, 75.10.Jm 



I. INTRODUCTION 

The discovery of high T c superconductors 20 years agoi 
led to an intensive search for new mechanisms, which 
could explain the observed superconductivity in GuOi 
planes. It became clear very soon, that two ingredients 
are needed, holes in low concentration which move in a 
2D antiferromagnetic background. The Hubbard model 2 
and the t — J model- have been discussed in this context. 
If the hopping of holes is frozen in, these models reduce 
to the Heisenberg model with Hamiltonian 



H = S(x)S(y) 



(1.1) 



where Si(x), i = 1, 2, 3 denote spin 1/2 operators at site 
x. In the absence of an external field (jl.ip conserves the 
total spin 



S = £S(*) 



(1.2) 



and the ground state is known to be a singlet state (S — 
0) with momentum p = (0,0). The first excited state 
\l,q) is a triplet (S = I, q = 0, ±1) with momentum 
p = (tt,tt). The transition between these two states is 
mediated by the staggered spin operator 



S * q = H(-l) x S,(aO 9 = 0,±1 



(1.3) 



So{x) = t^sOk), S±(x) = -^j=(ai{x) ±ia 2 (x) 



and the matrix element 



m 



N 



(U\s* q \o) 



(1.4) 



(1.5) 



can be considered as an order parameter. 

The ground state properties of the 2D Heisenberg 
model (|1.1|) have been investigated with various meth- 
ods. The variational RV B-stdX<& ("Resonating Valence 



Bond") starts from singlet states on pairs of sites, which 
cover the whole 2D lattice. By construction, these states 
lead to eigenstates of the total spin S with S = 0. How- 
ever, the manifold of these states which can be con- 
structed is nonorthogonal and overcomplete. 

Numerical methods - e.g. the Lanczos algorithm - arc 
limited to small clusters N < 6 x 6 = 36 due to storage 
problems. The computation on the largest cluster 6x6 
has been performed by Schulz and Ziman- 15 years ago. 

In spite of the great improvements achieved during this 
time, it is not possible so far to repeat the calculation of 
Schulz and Ziman for the next interesting cluster 8 x 
8. This is only possible with other techniques, as the 
quantum Monte Carlo [cf. e.g. (@)]. 

Our approach to the low energy properties of the two- 
dimensional antiferromagnetic Heisenberg model starts 
from a decomposition of the lattice into plaquettes as 
depicted in Fig. [TJ Quantum numbers, energies and cou- 
plings for plaquette states are discussed in Section [Til 
In particular we find, that the lowest energy plaquette 
states with total spin S = (singlet), S — 1 (triplet), 
S = 2 (quintuplet) appear to be most important for the 
construction of the ground states on larger clusters. This 
is demonstrated in Section IlIII where we compose 4 pla- 
quettes to a 4 x 4 cluster with open boundary conditions 
(cf. Fig. [2]). We deduce interaction matrices from the 
couplings between neighbouring plaquettes for the 4x4 
system in the sectors with total spin S = 0, 1, 2. The di- 
agonalization of these interaction matrices yield predic- 
tions for energies and transition matrix elements, which 
can be compared with Lanczos results on a 4 x 4 sys- 
tem. We find agreement within 1 — 6%, depending on 
the quantity under consideration. 

In a next step, described in Section IIV1 we generalize 
this approach to larger "n-clusters" of size 2™ x 2™, n = 
2, 3,4,.. .. E.g. the n — 3 cluster of size 8 x 8 is composed 
from four n = 2 clusters (4x4). Each of these clusters 
can be occupied with a cluster ground state with total 
spin S (S = 0, 1, 2). The interaction matrix for the n + 1- 
clusters have the same structure as in the step before - i.e. 



2 



for n-clusters. The n-dependence can be absorbed in a 
renormalization of couplings and gaps, which is discussed 
in Sections |V] and IVII The numerical evaluation of the 
renormalization group equations is discussed in Section 
IVIII Section IVIIII is devoted to the treatment of the 
staggered magnetization in our approach. 

Let us finally mention, that the idea to describe the 
ground state properties of the antiferromagnetic Heisen- 
berg model in the framework of a renormalization group 
approach is not new. One of the first attempts in this 
direction has been presented already in 1992 by Lin and 
Campbell.I^ They started from L x L clusters with L 
odd (L = 3, 5) and computed the ground state (which 
has total spin 1/2) and its interaction with sites on a 
ring. In this way, they were able to make a prediction for 
the staggered magnetization on a 7 x 7 lattice. 

There have been many investigations of ordered an- 
tiferromagnets, which start from the observation that 
the dominant fluctuations are controlled by the quantum 
nonlinear er-model with imaginary time^ ' 10 ^ 1 ' 12 

The various renormalization group approaches differ 
in the clusters used and in the truncation of the Hilbert 
space, which is needed to make the evaluation feasible. 
This is discussed in Section ILXl 



II. PLAQUETTE STATES: QUANTUM 
NUMBERS, ENERGIES AND COUPLINGS 

Our approach to the 2D antiferromagnetic Heisenberg 
model starts with a decomposition of the 2D lattice into 
plaquettes as depicted in Fig. [T] Each plaquette can 
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FIG. 1: 4-plaquette system with 2"x2"=4x4 sites - here 
n = 2; the single plaquette on the right shows the enumeration 
of plaquette sites. 



sentations of the SU(2). In this way we get two singlets 
|0), |0) with energies Eq — —2, Eq — 0, three triplets 
\Aq), Bq), \Cq) with energies E A = -1, E B = E c = 
and one quintuplet \Q r ) with energy E% = 1. 

The ground state |0) is given by the following spin 
configuration on the plaquette 




The three triplet states are obtained 

1 



- 2 



+ 



(2.1) 



\i,q) = SW|0) 



" ^(o|si l '4 l) |o) 1/2 



A,B,C (2.2) 



by application of the plaquette spin operators 

= I> (i) C«OS,(aO ? = 0,±1 i = A,B,C 



x£P 



(2.3) 



on the ground state (|2.1[) . The signs x^ l \ x ) 

<W{x) = (+ + —), 



x (A H x ) = (+-+-), 

X {C) (x) = (+--+) 



(2.4) 



define the magnetic order of triplet states on the plaque- 
tte: 

\A q ) is antiferromagnetic in the sense, that it changes 
sign running around the plaquette 

\B q ) and \C q ) are collinear antiferromagnetic in hori- 
zontal and vertical direction, respectively. 

The tensor states \Q r ), r = —2,-1,0,1,2 are totally 
symmetric with respect to the four sites of the plaquette. 

The simplest variational ansatz on the 2D lattice would 
start from a product state where all plaquettes are oc- 
cupied with singlets. Such an ansatz would lead to an 
energy per plaquette E (P) = —2 which is just 75% of 
the "exact" value 



carry 16 eigenstates. Their quantum numbers and ener- 
gies are listed in Table U 
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TABLE I: Quantum numbers of the 16 plaquette states. 



The spin quantum numbers S(P) in the first row result 
from a decomposition of the 4 spins into irreducible repre- 



E Q = 4 • (-0.668) = -2.674 (2.5) 

as it follows for the thermodynamical limit from a finite- 
size scaling analysis [HuseW (1988)]. Recent quan- 
tum Monte Carlo calculations [Sandvik 15 (1997), Loew-i£ 
(2007)] improve the ground state energy to eo = 
0.669437(5). 

Therefore the interaction between neighbouring pla- 
quettes Pi P r - as depicted in Fig. [2]- has to account for 
25% of the ground state energy. If the two neighbour- 
ing plaquettes carry spins \Si,mi), \S r ,m r ) S = 0,1,2, 
m = —S, . . . , S the spin interaction term S(x)S(y) at two 
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FIG. 2: Interaction of neighbouring plaquettes. 



neighbouring sites x G Pi, y G P r induces a change in the 
spin quantum numbers 



\Si,mt) S(x) \SiM) 



\S r 



| S J, - TYljy, i 



(2.6) 



[ si 


1 


* 


W 


<1 


mi J 




The transition matrix element 
{S' u m'i\S q {x)\S u mi) = 



can be evaluated by means of the Wigner-Eckart 
Theorem.— All these matrix elements can be expressed 
in terms of a Clebsch-Gordan coefficient and one reduced 
matrix element M(S[, x, Si). The latter only depends on 
the initial and final plaquette spins Si and and the 
triplet operator S q (x) at site x. The phase v q (v + = — 1, 
vq = V- = 1) results from the transformation proper- 
ties of the spin operator S q (x) under the group SU(2). 
The interaction between neighbouring plaquettes then 
depends on the product of two Clebsch-Gordan coeffi- 
cients 



v„v„ 





1 
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S r \ 


rn'i 


<1 


mi J 


V m 'r 
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m r J 



M(Sl,Si;S' r ,S r ) 



(2.8) 



and the product of reduced matrix elements 
M(S' h Si;S' r ,S r ) = ^2M(Sl,x,Si)M(S' r ,y,S r ) 

(2.9) 

summed over the two neighbouring sites (x,y), which 
connect the left and right plaquette as shown in Fig. [5] 
The Clebsch-Gordan coefficients lead to selection rules 



S'i — Si ± 1, Si 
S' t = 1 
mj = mi + q . 



for Si > 1 
for Si = 



(2.10) 



The explicit calculation of the transition matrix ele- 
ments 1 yields the weights M(S[, Si; S' r , S r ) for the fol- 
lowing cases: 



1 Note, that we always and without loss of correctness define the 
reduced matrix elements such that M(S',, Si; S' r , S r ) only depend 
on SI > Si, S' r >S r . 



1. Creation of triplet pairs from singlet pairs 

^ A q A- q (2.11) 



Note, that the initial and final states have total spin 
0. According to (f2T5)l and (|23|) the weight for this 
process turns out to be 



E 



where 
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9. 






-q I 



M (1) (10,10) 



M (1) (10,10) = = -i 



(2.12) 



(2.13) 



2. Hopping of an isolated triplet on a singlet back- 
ground 



A q <-> 



A q g = ±1,0 



(2.14) 



Here, the initial and final states are triplets on dif- 
ferent plaquettes (P r and Pi). The corresponding 
weight is given again by (|2.13p . 



3. The triplet-triplet process 

A s A t g = ±1,0 
introduces a new weight: 



A S + (] Af — g 



(2.15) 



i 

s + q 




where 



M (1) (ll,ll) 



1 1 

-q t 



M (1) (ll,ll). 



4. The process 



Q r 



<— ► A r — g Ag 



(2.16) 



(2.17) 



(2.18) 



starts from an initial state with a spin 2 plaquette 
Q r and a spin plaquette (singlet). In the final 
state we have two triplets A r - q A q coupled with 
an appropriate Clebsch-Gordan coefficient to form 
again a spin 2 state. The coupling for (|2.18[) turns 
out to be 




q ) 



M(21,10) 



where 



M (1) (21,10) = 1 



2\/3 



(2.19) 
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5. The process 



A r _ g 



describes the transition from a spin 2 - spin 1 to a 
spin 1 spin plaquette. It is accompanied with a 
weight 



{-) q 




1 



1 ° M ( V,10) 
-q I 



Truncating the states with subdominant transitions 
anticipates antiferromagnetic order of the system for the 
(2.20) renormalization procedure. This is somewhat in analogy 
to treatments by spin-wave theory or mappings to non- 
linear sigma models, see e.g. (fill ) and (jtj). In these ap- 
proaches long-range antiferromagnetic order is assumed 
from the beginning and fluctuations around this is built 
in subsequently. In our approach however, the system 
may or may not develop long-range order. This is deter- 
mined by the renormalization group flow. 



(2.21) 

6. The exchange of a spin 2 - spin 1 plaquette 

Q r Ap <— ► A r — q Qp+q (2.22) 

carries a weight 



2 




2 


1 


i 


r 


q r-q J 


\p + q 


q 





M (1) (21,21). 



III. THE FOUR PLAQUETTE SYSTEM 

We are going to construct in this Section the ground 
states of the four plaquette system with 4x4 sites de- 
picted in Fig. [T] It turns out that the ground states 
are symmetric under rotation of the four plaquettes. We 
therefore start from rotationally symmetric basis states, 
which are eigenstates of the total spin squared S 2 and its 
(2.23) 3-component 



where 



M (1) (21,21) = -i 



(2.24) 



(3.1) 



7. If the spin 2 - spin 1 plaquettes do not change their 
position 



Qr Ap <— > 

the corresponding weight is 



r+q Ap—q 



(2.25) 



2 

r + q 




1 1 » M (1) (22,ll). 



-q p 



where 



M (1) (22,ll) = ^. 



(2.26) 



(2.27) 



(|2 . 1 3|) - (|2 . 26[) is a complete list of processes where only 
singlets |0), A-triplets \A q ) and quintuplets \Q r ) are in- 
volved. In our opinion they are most important for the 
low energy behaviour for the following reasons: 

• |0) and \A q ) have the lowest energies according to 
Table H 

• the staggered spin operator [eq. (|2.3[) for i = A] on 
the plaquette induces the dominant transitions 



?(*) 



10} - \Aq) 



\Qr) 



As a consequence, the weights (|2.13l) . (|2.19p and 
(I2.24[) which are built up from these transitions are 
dominant as well. 



A. The singlet sector 

In Table |TT] we list 7 singlet states which can be con- 
structed on the four plaquette system with singlets (0), 
triplets (A q ) and at most one spin 2 (Q r ) plaquette. 



1 1,0) 
|2,0) 








4,°>=l £<->"+• (^i: 





A q 



+ 



rot 

A a 



|5',0> = E(-) p+9 < 




\ 


i6,o) = -^=Y(-r +q 1 2 


<7 P J 



A q Q. 



, + I | + rot 

/ V A v 



Q — r Ap 



+ rot 



TABLE II: Basis states of the 4 plaquette system in the singlet 
sector. 
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Starting from the state |1,0), where the four plaque- 
ttes are occupied with plaquette singlets |0), the creation 
process (|2.11|) generates the state |2,0) with neighbour- 
ing triplets A q A- q coupled to a total spin 0. The hopping 
process (|2.14l) leads from |2, 0) to |3,0). Further appli- 
cation of (12.111) and (|2.14|) generates the states 1 3, 0) , , 
|4,0), |5',0). The first four states |i,0), % = 1,2,3,4 are 
orthonormal. This is not the case for |5',0), which is 
orthonormalized by 



[5,0) = A%(|5',0)-<4,0|5',0)|4,0)), (3.2) 



where 



and 



(4,0|5',0) 



(3.3) 



2. The triplet-triplet and quintuplet-quintuplet tran- 
sition matrix elements, which define the weights 



/3« 



MW(11;11) _ 3 
tf) ~~4 

M«(22;ll) _ 3\/3 
TJX) ~ 4~ 



(3.11) 
(3.12) 



only contribute to the diagonal matrix elements. 
They also depend on the two scaled gaps 



P = P 



i) _ 1 _ o 



(3.13) 



(3.14) 



4^5 



(3.4) 



The states |6, 0) and |7, 0) contain one spin 2 plaquette 
Q r coupled together with two spin 1 plaquettes A q A p 
to form a state with total spin 0. The states \i, 0) i = 
1, 2, .., 7 are orthonormal and complete in the sense, that 
no further rotational symmetric state can be constructed 
from plaquette singlets, triplets and one quintuplet. The 
Hamiltonian restricted to the Hilbert space of these seven 
singlet states can be written as 



(3.5) 



The first term is just the energy of the state |1,0). We 
have scaled out the singlet-triplet coupling (|2 . 13f) . The 
remaining "interaction matrix" Aq is listed in Appendix 
IAT1 

The following remarkable features can be observed in 
the interaction matrix A 



(2). 



1. The nondiagonal matrix elements are nonnegative 
and fixed by the weights 



= M (1) (10;10) 

M (1) (21;10) 
M (1) (10;10) 



1 

3 

V3 
2 



(3.6) 
(3.7) 



They are induced by the singlet-triplet (1,0) and 
triplet-quintuplet (2, 1) transition matrix elements 
according to (|2.9[) . 

Therefore, the Perron-Frobenius theorem holds, 
which states that the eigenvectors \cr^) with 
largest eigenvalue a^: 



<rW\*W) 



cr (2) = 2.20917 
have nonnegative components: 

M|ct (2) ) > i= 1,...,7. 



(3.8) 
(3.9) 

(3.10) 



3. The ground state energy of the Hamiltonian 
in the restricted Hilbert space of the singlet states 
|i,0), j = 1, ...,7 is given by 



a * (j = -8.7236 (3.15) 



which deviates from the "exact value for the 4x4 
system with open b.c. 



by 5.0%. 



E^ L) = -9.1892 



B. The triplet sector 



(3.16) 



The rotational symmetric eigenstates of the 4 plaque- 
tte system in the sector with total spin S = 1 are listed 
in Table [ml 

Starting from |1,+), we generate the other states 
|3', +), |2, +), |4, +) and |5, +) by means of the processes 
(|2~TT1) (|2~T41) and (|2~TSD . The states |1,+), |2,+), |4,+), 
|5,+) arc orthonormal. The state |3',+) is not yet or- 
thonormal with respect to |2, +), which is achieved by: 



IV 



N 3 (\3',+) - (2,+|3',+)|2,+)) , (3.17) 



with 



and 



<2, + |3', + > = ^ 



V3_ 

4\/5 



(3.18) 



(3.19) 



The states \k,+), k = 6,7,8,9 contain one spin- 
2 and three spin-1 plaquettes, which are coupled to- 
gether with appropriate Clebsch-Gordan coefficients 
(CL'r _ p _g) eigenstates with total spin 1. There exist 
two further states (j = 1, J = 1,2) that do not couple to 
the considered ones. 



6 



n,+>= 4 \ A+ 

v 7 ! I v 



rot 

A+ A 



+ rot 



i3',+> = E(-) 9 { 












2 1 >< 

i-9 9/ 



|5,+> = 



1 



3l-9 j4<j 


2 1 
1-99 



+ | " | + rot 



,, 







A 



+ rot 



\k 
with 



J.J 

-r,—p, —q 

J 



3 

■p-q 

J 1 
p — g — r + p + g 



1 1 

-p -g 



X 

2 J 
r + 1 — r 





k 


6 


7 


8 9 


and: 


j 





2 


2 2 




J 


1 


1 


2 3 



TABLE III: Basis states of the 4 plaquette system in the 
triplet sector. 



The Hamiltonian in the restricted Hilbert space of the 
states \k,+) k = 1, . . . , 9 reads 



(3.20) 



The first term is the energy of the lowest state 
Again we have scaled out the singlet-triplet coupling 

(|3.6|) . such that the interaction matrix depends on 
the two scaled gaps ([3~T3]) . (|3~T4|) and the three 
scaled couplings 7 « (|577jl . /?W (13~TD) and e« (l3~T2l . 
Diagonalizing the interaction matrix: 

Af } |T (2) > - r^lr^) (3.21) 
yields the largest eigenvalue 

r (2) = 3.41009 (3.22) 
which corresponds to a ground state energy 

El 2 ' = -7- — = -8.18405. (3.23) 

The latter deviates from the Lanczos result on a 4 x 4 
system with open b.c. 



E[ L) = -8.6869. 



(3.24) 



by 6.3%. 



C. The spin 2 sector 

The 14 basis states \l, 2+) for the 4-plaquette system 
(Fig. [T|) in the spin 2 sector are listed in Table HVl 

The states \l, 2+), / = 1, 2, 3, 5 are orthonormal, which 
is not the case for |4',2+) with respect to |3,2+). We 
therefore introduce 



|4,2+) = A^ 4 (|4',2+) - (3,2 + |4',2+>|3,2+>) 



with 



(3,2 + ]4',2+) = 
N 4 = 



4^/2 



(3.25) 

(3.26) 
(3.27) 



The states \6+j, 2+), \9+j, 2+) and |12+j, 2+) contain 
one spin 2 and two spin 1 plaquettes - the latter are 
coupled together to a state with spin j, which then forms 
with the Q r plaquette a state with total spin 2. 

In the restricted Hilbert space of the states \l, 2+), I = 
1, . . . , 14 the Hamiltonian can be written as 

Bf = 2(4 1} Af . (3.28) 

Again the first term corresponds to the plaquette energies 
of the state |1,2+) (and |2, 2+)). The interaction matrix 

A< 2) depends on the scaled gaps (|3~T3j) . «W 
and the three scaled couplings (|3.7p . /3W p. lip and 
e' 1 ) (|3.12p . as can be seen in Appendix I A 31 
Diagonalizing the interaction matrix 

A^ 2) |C (2) ) = C (2) le (2) ) (3.29) 
yields the largest eigenvalue 

£ (2) = 3.48987 (3.30) 
which corresponds to a ground state energy 



Ei 2) = -6 - ^— = -7.36767 . 



(2) 



(3.31) 



The latter deviates from the Lanczos result on a 4 x 4 
system with open b.c. 



E. 



(L) 



-7.7909. 



(3.32) 



by 8.0%. 



IV. 



THE RENORMALIZATION GROUP 
PROCEDURE 



In the previous Section we have explained how to con- 
struct the ground state on a (n = 2) cluster of size 4x4 
from four interacting n = 1 plaquettes (2 x 2). We only 
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|5,2+> = 
|6+j,2+) = 
\9 + j,2+) -- 
|12+i, 2+) : 



1 
1 

71 



2 2+ o 
o o 



+ rot 



74 



r<3 

^p,q — 



A — q A — p 
Qr + 2 

A- p 

Qr+2 A- q 

2 3 
2 + r -r 



rot 



+ rot 



+ rot 

1 1 

-q -p 



j = 0,1,2; r=p + q 



TABLE IV: Basis states of the 4 plaquette system in the quin- 
tuplet sector. 



took into account plaquette states with total spin 0, 1 
and 2. This procedure will now be extended to compute 
the ground states on (n+ 1) clusters (2™ +1 x 2" +1 ) with 
total spin 5 = 0, 1,2 from the corresponding quantities 
on n clusters. The ground states 

|5,m;n+l) S = 0, 1, 2, m = -5, . . . 5 (4.1) 

are supposed to carry alternating momenta 

p s = (0, 0) for 5 = 0, 2 
Ps = (n, 7r) for 5 = 1. 

(4.2) 

Our approach starts from the basic assumption, that the 
ground state (|4.1[) on (n + 1) clusters can be constructed 
from the ground state on n clusters 2 

\S',m';n) S = 5' + 1, 5', 5' - 1 5' > 1 
5=1,0 5' = 

(4.3) 

Under this assumption the Hamiltonians Hg l+lS> on 
the (n + 1) cluster (2" +1 x 2 ,l+1 ) can be written in an 
analogous form to J33]) for 5 = 0, (j3T2"0|) for 5 = 1 and 



2 Of course the full Hilbert space contains many more states on n- 
clusters — like higher spin states and states with momenta p ^ ps 
like the B and C triplets on the n = 1 plaquette (Table [IJ. 



(|3.28p for 5 = 2. To be definite, we introduce first the 
analogues of the basis states in Tables HU Mil II VI 



\i, 0; n H 


-1) 


i = l,. 


.7 


(4.4) 


\k, + ;n- 


-1) 


k = 1,. 


..9 


(4.5) 


l,2+;n- 


-1) 


1 = 1,.. 


.14 


(4.6) 



tn+1) 



on a (n + 1) block (2 n+1 x 2" +1 ). 

Then the analogues of the interaction matrices Aq 
are obtained from Appendix [Al by substituting the scaled 
energy gaps and couplings 







^0 


(4.7) 


a («) 


«<"> 






(4.8) 




a (") 




= M {n) 


(10; 10) 


(4.9) 


7 («) 


1 

~ aW 


•M (n) (21;10) 


(4.10) 




1 

~ aW 


•M (n) (ll;ll) 


(4.11) 


£ (n) 


1 


•M (n) (22;ll) 


(4.12) 



The latter can be related by 

M (5;, Si', S r , S r ) 



V^M(")(5 z ',x,5)M(")(5;,y,5 r ) 



to the reduced matrix elements 

{S[,m' l ,n\S q {x)\Si,mi,n) 

M^\S[,x,Si) 



(4.13) 



* 


1 




V m 5 


9 





(4.14) 

for the transition Si, mi — > 5j,mJ of the cluster spins. 

The renormalization group procedure only affects the 
reduced matrix elements - i.e. the couplings (|4.9p - (|4.12p 
and the scaled gaps (|4. T|) and (|4. 8|) which enter as pa- 
rameters in the analogues for the interaction matrices in 
Appendix [X] 

A^ l+1) (pW , K («) , 7 (») , /3(«) , e («) ) 5 = 0,1,2 

(4.15) 

on an (n + 1) block 2 n+1 . The largest eigenvalues of the 
interaction matrices 



.(n+l) | „(«+!) \ 



c7 („+i),(„+i) ) 



(4.16) 



A («+D| r (r l+ l) ) = T (n+l)| T (n+l)j (417) 
_ £(n+l)|£(n+l)\ (4.18) 



A (n+l)|^(n+l) 
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yield for the ground state energies 





= AE^ ] - 


ha (") CT (n+l) 


(4.19) 


E (n+1) 


= 3E { n) - 


- P< rl) + a ("M" +1 > 


(4.20) 


E (n+1) 


= 2(4 n) 


+ E[ n) )+a (n) ^ n+1) . 


(4.21) 



V. THE RENORMALIZATION OF THE SPIN 
MATRIX ELEMENTS 

Our starting point is the group of spin matrix elements 
on a (n + 1) cluster: 

(T (n+ V\S + (x)\a {n+1) ) = M (n+1 \l,x,0) (5.1) 
(£(«+i)|S + ( x )| T (»+i)) = M(™+ 1 )(2,a;,l) (5.2) 



(t^\S (x)\t^) 
(& +1) \S Q (x)\& +1) ) 



1 1 

1 

1 2 
2 



M (n+1) (1,3,1) 
(5.3) 

M(" +1 )(2,a;,2) 
(5.4) 



These matrix elements are expressed in terms of the cor- 
responding reduced matrix elements by means of the 
Wigner Eckart Theorem [cf. f27T|)]. 

The eigenstates of the interaction matrices [cf. (I4.16[) - 
(I4.18|) ] are represented in terms of the basis states [ (I4.4[) - 



r (n+l)\ = 



J2^r +1} \iAn+l) (5.5) 

i=l 

|r ( " +1) ) = ^rf +1 ^,l,n+l) (5.6) 



fc=i 

14 



= X;e I (B+1) |i,2,n+l) (5.7) 



1=1 



which leads to the following set of recursion formulas for 
the reduced matrix elements 

Af( n+1 )(l,a;,0) = Z( n+1 >(l,0)M( n > (l,a:,0) 

+/(»+!) (2, 1)MW (2, s,l) 
M(" +1 )(2,a;,l) = G (n+1) (1, 0)Af (n) (l, x, 0) 

+G (n+1) (2,l)M (n) (2,x,l) 
M(" +1 )(l,a;,l) = F(" +1 '(1,1)MW(1,i,1) 

+i^(«+i)(2,2)M^(2,a;,2) 
M(" +1 )(2,a;,2) = F ? (n+1) (l, 1)M ( ")(1, a:, 1) 

+F 5 (n+1) (2,2)M(")(2,x,2). 

(5.8) 



The coefficients depend in a bilinear form on the compo- 
nents of the eigenvectors 

J(" +I )(a,6) = 5>*" +1)j MM)^ n+1) (5-9) 

fc,Z 

G(" +1 )(a,6) = ^^ +1) G ; , fc (a,6)rf +1) (5.10) 

F^foo) = £(r(" +1) ) 2 F T , fc ( a ,a) (5.11) 
fc 

Ff+^a) = ]T(^ +1) ) F^(a,a). (5.12) 
z 

Nonzero elements of the total spin combinations (a,b) — 
(2,1) and (a, a) = (2,2) are marked with boxes in the 
Tables of Appendix IBl 

Note that the renormalization - i.e. the n-dependence 
of the coefficients (|5.9p - (|5.12p - only enters via the eigen- 
vector components. The "contractions" 7,^(1,0), etc. 
are independent of n and solely determined by the spin 
matrix elements between the basis states [ (|4.4I) - (|4.6[) ]. 
Therefore, they have to be calculated once and are listed 
in Appendix [Bj 

We can check the quality of the recursion formulas 
[M(" +1 )(l,a;,0),M(™ +1 )(2,a;,l) - cf. ([STS]) ] in the first 
step (n = 1): 



7< 2 )(1,0) = 


-0.64796 


(5.13) 


/( 2 )(2,1) = 


+0.03387 


(5.14) 


G (2) (1,0) = 


-0.66177 


(5.15) 


G (2) (2,l) = 


+0.02062 


(5.16) 



if we compute the transition matrix elements for the stag- 
gered spin operator S* (P) on a 2 x 2 plaquette 



M (1) (1,P,0) = p = -1.63299 (5.17) 

V6 

M (1) (2,P,1) = —= = 1.41421 (5.18) 
2v 2 

Af (2) (l,P,0) = +0.64796 ■ —= + 0.03387 • V2 

V6 

= 1.10601 (5.19) 

M (2) (2,P, 1) = +0.66177- —= + 0.02062- V2 

V6 

= 1.10983 (5.20) 

and compare it with the Lanczos result on a 4 x 4 system 
with open b.c. 

M (i) (l,P,0) = 1.0857 (5.21) 
M (i) (2,P, 1) = 1.1826 (5.22) 



(j5TT§|) and (jST2"0|) deviate from (IOT|) and by 
+1.9% and —6.2%, respectively. 
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VI. RECURSION FORMULAS FOR THE 
SCALED COUPLINGS AND GAPS 



The relevant couplings between neighbouring n-blocks 
as depicted in Fig. [3]- can be expressed via eqn. (|4.13|) 




FIG. 3: Couplings between neighbouring n-blocks - here 
shown for n = 2. 



in terms of the reduced matrix elements defined through 
(|4.14j) . In (|4.13|) we have to sum over the 2" bonds 
(x, y) , which connect the left and right block. The nearest 
neighbour interaction S(x)S(y) on these bonds changes 
the total spin on the left and right block 



Si 



s; 



S ., 



in the same way as we discussed in (12.61) for the plaquette 
interaction (n = 1) as depicted in Fig. [2] Note that the 
definitions (|4.13p and (|2.9p are identical for all blocks of 
sizes 2™ x 2™, provided we perform the summation over 
the 2™ connecting bonds correctly. 

If we admit only blocks with total spin 0, 1, 2 we get 
from (|4.13|) and the recursion formulas (|5.8|) the renor- 
malization of the couplings (|4.9l) - (14.12j) : 



q("+i) 

/«+!) £: 1 

1 2oW 



(V" +1 )(l,0)+7 (rl) / (n+1) (2,l)) 

(6.1) 

^ / (™+i)( lj 0)+ 7 (")/(»+ 1 )(2,l)) 

•(G ( " +1) (l,0) + 7 (n) G ( " +1) (2,l)) 

(6.2) 

i?(™+l)(l ) l)2. /3 (n) + 

2Fj n+1 \l, l)F^ n+1 '> (2, 2) • e<"> (6.3) 
F^ +1 )(l,l)if l+1) (l,l) -(3 (n) + 
>i n+1 )(l,l)Ff +1) (2,2) + 
i^ n+1) (l,l)F^ n+1 )(2,2)l -e {n) (6.4) 



In addition to the scaled couplings 7 < - n \ (3^ and the 
interaction matrices (|4.15p depend on the scaled energy 
differences (I4T71) . (gTSJ). From (|4TT9"|) - (|4T2"T]) we get the 



fin+l) 



2a(") 

q("+ x ) 
2a(") 



recursion formulas 

p(»+l) a (n+l) 
q(n) 



T (n+1) _ CT (n+l) 

'-0 ^ ,(, ' : " 



/0 



(6.6) 

Here, [J] denotes the positions where a plaquette- 
plaquette interaction (of general coupling strength J) 
would have to be implemented. Throughout this work, 
however, we will keep J = 1 and discuss the interesting 
case of a phase transition in a two-dimensional model of 
interacting plaquettes in a separate paper (|24j). 

In summary, we see, that each step n — > n + 1 in the 
renormalization procedure demands the diagonalization 



of the three interaction matrices A 



(n+l) 



S = 0,1,2. The 



eigenstates |cr(™ +1 )), |r( Tl+1 )), |£( Tl+1 )) with largest eigen- 
values a(' n+1 \ r( n+1 \ determine the right-hand 
sides of the recursion formulas (|6.ip - (|6.6p . 



VII. NUMERICAL EVALUATION OF THE 
RENORMALIZATION GROUP FLOW 

In this section we present numerical results for the evo- 
lution of couplings and scaled gaps with n, which defines 
the block size 2" x 2™. We start from the states in Tables 
HTl IIIIl ITVl for the singlet, triplet and spin-2 sectors. The 
dimensions ds of the interaction matrices As S = 0, 1, 2 
increases with S: 



do = 7, 



di = 9, 



14 



(7.1) 



since the number of possibilities to construct 4 plaque- 
tte states with singlet, triplet and at most one spin 2 
plaquette increases with 5*. 
In Fig. HJ the ratio 



2a(") 



a < 0.525 for n > 5 



(7.2) 



is shown; it converges to a constant value slightly above 
1/2- 

In Fig. [5j we present the evolution of the couplings 
([621), eW (JEU) 

(n) 



7* = 1.0847 (7.3) 
(7.4) 
e^-' -> (7.5) 

This result has to be interpreted that the couplings for 



7 
fin) 



• the "nondiagonal" transitions (|4.10p . (|4.9p with 
spin exchange 

Q A «-► 4 , ^O^OA 

are both relevant in the vicinity of the fixed point. 
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2 4 6 8 10 12 14 16 18 20 




10 12 14 16 18 20 



FIG. 4: n-dependence of the coupling ratio " 2 (n) for dimen- 
sions do = 7, di — 9, d2 = 14. 



FIG. 6: Evolution of the scaled energy differences p^ n \ K^ n ' 
and of the largest eigenvalues cr (n+1) , r (n+1) , ^ n+1 >. 




4 6 S 10 12 14 16 18 20 



FIG. 5: Evolution of the couplings 7 (n+1) , /3 (n+1) , e (n+1) for 
dimensions do = 7, di = 9, d2 = 14. 



• the "diagonal" transitions (|4.1ip . (|4.12p with no 
change in the plaquette spins, however, 



A A <-» A A , Q A <-» Q A 



are irrelevant for n — ► oo. 



Note that the ratio 7( n+1 ) (|6.2[) only shows a slight 
variation between 1.086 and 1.151 which implies that 
the nondiagonal elements in the interaction matrix Ag, 
S = 0, 1,2 are almost constant. The diagonal elements 
depend on the scaled energy differences p^ (|6.5p and 
«W (|6.6p which increase with the system size 2™ x 2™ , as 
is shown in the lower part of Fig. [6] From the upper part 
we see that the largest eigenvalues a^ n+1 \ T^ n+1 \ 
of the interaction matrices increase with n. Indeed the 
essential mechanism of the renormalization group con- 
sists in a feedback between the increase of the (negative 
valued) quantities and the largest eigenvalues 

cr ( n + 1 ) ; T (™+ 1 ) j £(«+!). 



For large n 



p (n+l) _ 




= -0.046 


K (n+1) _ 


■+ 


= -0.368 


a (n+D _ 


-> <7* 


= 6.53 


r («+l) _ 


-> r* 


= 6.53 






= 6.15 



(7.6) 
(7.7) 
(7.8) 
(7.9) 
(7.10) 



the fixed point values p* , k* are close to zero, whereas 
a*, t* , £* approach each other. The deviations are a 
consequence of the reduced dimensions do = 7, d\ = 9, 
d,2 = 14 of the Hilbert spaces for the interaction matri- 
ces. We expect that these deviations will be lowered, if 
we enlarge ds, S = 0,1,2 systematically such that the 
energy differences 



E 



(«) 



E. 



(n) 



E 



i 

(n) 



2E 



(n) 


(n) 



^—nvi 
A—nv 2 



vanish in the thermodynamical limit n 
nents 



(7.11) 
(7.12) 

oo. The expo- 



log(l + x) 
log 4 



V2 



log(l + y) 
log 4 



(7.13) 



can be determined from the first derivative of the eigen- 
values cr, r, £ with respect to p and k, respectively: 



d{r 



din- 



dp 

d{a + ^-2r) 



dp 



8{t — a) dn 
Ok dp 



(7.14) 



d((T + £- 2t) 
8k 

| d(a + £-2T) 
dp 



dp 
dn 



(7.15) 



The partial derivatives of the eigenvalues a, r, £ with 
respect to the parameters p and k, which enter linearly 
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in the diagonals of the interaction matrices Ag(p, re) (cf. 
Appendb JA"]) can be computed from the matrix elements 
of S = 0,1,2 between the eigenstates \a), |r), 

10 (cf. Appendix E}. 

% = {a\ d -^-\a)=A{l-al)-2{al + al) 

Tp = H^l->=4(l-r 1 2 )-2(r| + r 3 2 +r|+r|) 

| = ( ^) = 2.(i-^-^-e3 2 ) 



r (n+l)\ 



5.5[) which enters in the definition of the stag- 



,0A Ol 



(J = (J fi 



\ -i 2 2 2 

r) = 1 - T - T, - To 



da 

dn x " 1 dre 
dr 8Ai 

dn {Tl dn 
1 = ^)^1-^1-^1 

(7.16) 

Remember, oi (i = 1, . . . , 7), Tk (k = 1, . . . , 9), (Z = 
1, . . . , 14) denote the components of the eigenvectors \a), 
\t), as they follow from the diagonalization of the 
interaction matrices As(p, re) for p — ■> 0, re — > 0. 




70 72 74 /6 75 20 



FIG. 7: Quantities a;,y [cf. eqns. (|7.14[1 , (|7.15|) ] for dimen- 
sions: do = 7, di = 9, d,2 = 14. 

In Fig. [7] we have plotted the quantities x and y vs. n 
which enter in (fTTH)) and (f7TT5|) . 

The quantity a^™) appears to be rather stable around 
-0.75 for n > 3 and yields a value of v\ = 1 for the crit- 
ical exponent appearing in (|7.1ip . The vanishing of the 
singlet-triplet gap (|7.1ip in the thermodynamical limit 
has been suggested by Tang and Hirsch (|18l ) from a finite- 
size analysis of the ground state energies. On the other 
hand is not yet stable with respect to n. Again the 
reason might be that our truncation of the dimensions 
ds for the interaction matrices As S = 0,1, 2, still needs 
to be improved. 



VIII. THE STAGGERED MAGNETIZATION 

The real space renormalization group approach gen- 
erates a recursion formula for the singlet ground state 



gered magnetization 



(n+l)|y(«+l)y(n+l)| (n+l)\ 



Here 



^O+i) 



X 



(8.1) 



(8.2) 



i=i 

is the properly normalized staggered spin operator on a 
(n + l)-cluster, which can be decomposed into the corre- 
sponding quantities on the four neighbouring n-clusters 
Pj, j = 1,2,3,4, as shown in Fig. [1] This leads to the 
following recursion formula for the ratio 



R (n+1) 



( a (n+l) | E (n+l) s (n+l) ^(n+l)) 
(a(«)|£ W £V"V n) ) 

= E ^ +1) ^ +1) ^Al [n) ) (8-3) 



In the evaluation of (|8 . 3[) we can use the fact that the 
singlet basis states |i,0) (Table ITT1) are invariant under 
rotations of the four plaquettes Pi, P 2 , P3, P4. Therefore, 
we are left with the computation of the matrix elements 



(i', 0;n+ 1|£W(P.)£^ ( Pl) | if 0; n + 1), j = 1, 2, 3 

(8.4) 

which proceeds in the following steps: 

a) The application of staggered spin operators 

(Pi) onto the singlet states \i,0;n+ 1) leads 
to triplet states, which are not rotational invariant. 
A convenient set of triplet basis states is defined in 
Appendix ICl [ (ICTfl - (|C7f ] together with the decom- 
position of 



sV" +1) (Pi)N,0;n+l) 



.5) 



into these basis states [ (|C10I) - (|C17|) ]. 
b) By rotation of the 4 plaquettes 

P1-P2 P 2 ^P 3 P 3 ^P 4 P 4 ^P (8.6) 



we obtain from the decomposition of (|8.5[) the cor- 
responding decompositions 

£+(P 2 )|i,0;n+l) £+(P 3 )|i,0;n + l) (8.7) 

into the triplet states |UT0|) - ([UT7)) . 

c) The scalar products (|8.4p turn out to be propor- 
tional to 

AfW(l, P, 0) 2 MW(1, P, 0)M (n) (2, P, 1) 
M {n) (2,P, if (8.8) 
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where 



M^(S+1,P,S) = ^(-) a Af< n >(S +!,*,£) 



xeP 



(8.9) 



are just given by the reduced matrix elements (|2.7|) 
summed over all sites of the n-cluster P. 



MW(1,P,0) 2 = (o-^lE^S^V^ 3 ) (8-10) 

can be identified with the staggered magnetization 
on the n-clustcr P, whereas the ratio 



M(")(2,P,1) _ (n) 



11) 



M(«)(1,P,0) 
is given by (|4.10p . 
d) The decomposition 

(*', 0; n + 1|S W (Pi)S^ (Pi)\i, 0; n + 1) = 

rJfAfWji, p o) 2 + r^; 2) M(")(2, p i) 2 

+r^Af P, 0)Af W (2, P, 1) (8.12) 

illustrates, that the n-dcpcndcncc - i.e. size depen- 
dence 2 n+1 x 2" +1 - only enters via the reduced ma- 
trix elements MW(1,P,0), Af( ,l )(2,P,l), whereas 
the 7x7 matrices 



r (j>o) r 0'.i) p0'.2) — i 7 j-ni 



(8.13) 



are independent of n. They are not affected by 
the renormalization group procedure and can be 
completely expressed in terms of scalar products 
formed from the triplet states (cf. Appendix [Cj) . 



This leads to an explicit expression (|C19|) - (|C22p of the 
7x7 matrix IV ^(7'™)), which enters into the recursion 
formula (|8.3[) . The numerical evaluation of (|8.3[) is pre- 
sented in Fig. [5J 

The ratio starts around 1/2 and rapidly increases to 
0.742 and remains constant for n > 4. A nonvanishing 
staggered magnetization would demand a limiting value 
R = 1 in the thermodynamical limit n —> 00. The de- 
viation from this value, we see in Fig. [5J has to be at- 
tributed to the truncation of the interaction matrices As, 
S = 0, 1, 2. Their dimensions ds are limited to 



(do,di,d 2 ) = (7,9,14) 



(8.14) 



since we allow only for one quintuplet plaquette on the 
four cluster compound. We expect that the extension of 
the interaction matrices to four cluster states with 2, 3, 4 
quintuplets will lead to an increase of the ratio R. For the 
moment, we can only compare the difference of including 
one quintuplet (|8 . 14[) to zero quintuplet contributions: 



(5,3,4) 
(7.9,14) 



2 4 6 8 10 12 14 16 18 20 



FIG. 8: Evaluation of the recursion formula (|8 . 3|) for dimen- 
sions: (d ,d 1 ,d 2 ) = (5,3,4), (7,9,14). 



(do,di,d 2 ) = (5,3,4) 



.15) 



The ratio R is substantially lower in the case (|8.15j) as 
can be seen from Fig. [5J 

In other words: Higher plaquette excitations are 
needed to generate plaquette-plaquette interactions 
which yield a nonvanishing staggered magnetization in 
the thermodynamical limit. 

It has been observed already by Bernu et al. (|19l ) that 
the ground states in the total spin S sectors collapse to 
the ground state in the thermodynamical limit. 



IX. DISCUSSION AND PERSPECTIVES 

If we compare our approach with previous renormal- 
ization group methods, we find on one hand the same 
goal - namely the derivation of a low energy effective 
Hamiltonian - but also crucial differences in the underly- 
ing assumptions: 

Most of the "older" approaches like that of Lepetit and 
Manousakis start with blocks with an odd number 
of sites. Here, the block ground state has spin 1/2 and 
the Wigner-Eckart Theorem allows already interactions 
between neighbouring blocks in the ground state. Ex- 
cited states - e.g. with total spin 3/2 - are assumed to 
contribute only to a renormalization of the coupling be- 
tween blocks in the ground state. Therefore, there is 
no renormalization of the energy difference between the 
ground state and excited states. In our opinion this is the 
reason, why these approaches do not allow for spin-spin 
correlations at large distances. The exact RG flow acts in 
an infinite-dimensional space of Hamiltonians resp. cou- 
plings. Even when starting with a model that is defined 
by very few couplings, the exact flow will carry the Hamil- 
tonian into rather complicated regimes: at each step of 
an RG procedure longer-ranged couplings are generated. 
In the past, in many applications of the RG concept the 
space of all Hamiltonians was truncated to a finite di- 
mensional one, i.e. only a few coupling constants were 
kept. For many universal properties, this approach was 
successful. 

Recent approaches like ours and that of Capponi et 
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al. 421D, Albuquerque et al. ((2J) based on the CORE 
method (Contractor Renormalization group) (23) start 
with plaquettes having a singlet ground state. They 
cannot interact, since the interaction is mediated by 
the spin operators in the Hamiltonian. Their expecta- 
tion values between total spin states vanish accord- 
ing to the Wigner-Eckart Theorem. Therefore excited 
states on the plaquettes are absolutely necessary to gen- 
erate interactions between the plaquettes [cf. the pro- 
cesses (I2^ . (l23il) . (l2^ . (l2^ . (l2~2^ . (l2~^2l . (l2"2D) of 
Section II]. For this reason we included triplet (|^4 g ), 
q = ±1,0) and quintuplet {\Q r }, r — ±2, ±1,0) exci- 
tations. Indeed it turned out that the states with one 
quintuplet excitation (i.e. |6,0) and |7,0) in Table HT1 and 
\k,+), k = 4,5,. .,9 in Table III] and |3,2+), |6 + j,2+), 
|9+j,2+), |12+j,2+) (j = 0,1,2) in Table EJ) improve 
the decrease of the singlet-triplet gap in the large-n limit 

(5, 3, 4) = -0.767 ; p* (7, 9, 14) = -0.183 (9.1) 



APPENDIX A: THE INTERACTION MATRICES 

Here we present the interaction matrices As, S = 
0, 1,2. In order to clarify the dependence on the scaled 
energy gaps p, k and coupling constants 7, j3, e [cf. eqs. 
(j4.7|) - (|4.12|) ]. it is convenient to consider the following 
block forms: 



1. Spin sector 



A = 



A(5,5) 
A T (5,6) 



A(5,6) 
A(6,6) 



(Al) 



A(5, 5) is fixed by the matrix elements of the first five 
states: 



The limiting value (|9.ip defines a measure for the "qual- 
ity" of the singlet-triplet gap generated with interaction 
matrices A5 of dimensions ds, S — 0, 1, 2 



(d ,d 1 ,d 2 ) = (5, 3, 4) ; (7, 9, 14) 



(9.2) 



Such an extension of the interaction matrices also leads 
to an improvement of the staggered magnetization, as 
discussed in Fig El 

In refs. (f2lh . ((221 ) the quintuplet excitations are miss- 
ing and it would be interesting to see how the singlet- 
triplet gap decreases in their renormalization process. 
Note however, that quintuplet excitations possess large 
couplings (|2.19p to triplet excitations, which do not die 
out in the renormalization process (n — > 00) (Fig. [S]); the 
corresponding energy differences (14. 7|) , (|4.8[) decrease as 
well (Fig. [|5]). The authors of ref. (|22j) intend to improve 
the CORE results by varying the compounds of plaque- 
ttes. In addition to the 4 plaquette compound of Fig. 
[3 they allow 2 and 3 plaquette compounds. We do not 
have this freedom, since our renormalization approach is 
restricted to the geometry of the 4 plaquette compound. 
The restriction to rotational symmetric states on the 4 
plaquette compound with singlets, triplets and quintu- 
plets enables us to follow the renormalization group flow 
for all couplings and gaps, which enter into the interac- 
tion matrices. 

We intend to improve our results in a first step by 
taking into account all rotational symmetric states on 
the 4 plaquette compound with uq = 2,3,4 quintu- 
plets. Larger interaction matrices demand for an efficient 
method to calculate matrix elements which form SU (2) 
invariants - similar to Racah coefficients (fl7T ) in Nuclear 
Physics. 
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Note that this block matrix only depends on p and p. 
The matrix elements between the states i = 1, . . . , 5 
and i = 6, 7 are contained in the 2x5 block matrix 

2,/| ^ ^ 
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which is proportional to 7. 

The matrix elements i = 6, 7, i' = 6, 7 form the third 
block matrix 
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2. Spin 1 sector 



Ai = 
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A(5, 5) is fixed by the matrix elements of the first five 
states: 
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A(5, 5) is fixed by the matrix elements of the first five 
states: 
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The matrix elements between the states k 



and k — 6, 7, 8, 9 are contained in the 4x5 block matrix 
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The matrix elements k 
the third block matrix: 



A(6,6) = (4p + «)-l 
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3. 


Spin 2 sector 
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A(6,9) = F + 7 2 C* 
A(6,12) = F + ~f 2 FCF 



(A20) 
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APPENDIX B: COEFFICIENTS (|5^) - (|532]) 



In this Appendix the coefficients (|5.9p - (|5.12p of Section 
|V]are given. Note that the dimensions ds — 7, 9, 14 have 



Note further, that those nonzero coefficients with spin 
2 plaquette contribution are marked with boxes. I.e., 
the reduced Hilbert space [(do, d\, d%) = (5, 3, 4)] for zero 
spin 2 plaquettes is directly illustrated. 
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APPENDIX C: RECURSION FORMULA FOR 
THE STAGGERED MAGNETIZATION 



(B3) 



As explained in Section IVIIII the derivation of 
follows the steps 



a) The application of the staggered spin operator 

£j"^(Pi) onto the singlet basis states \i, 0;n + 1) 
- listed in Table |TT] - generate triplet states, which 
are not rotational invariant: 
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It is convenient to express these states in a new 
basis which defines the occupation of the plaquettcs 
Pi P2 P3 P4 in cyclic order, e.g. 

e£°(Pi)|1,0) = -M(")(1,P, 0)|A+, 0,0,0) (C8) 
where 

-M<">(1,P,0) = (A f ,0,0,0|S+(P 1 )|l,0) (C9) 
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In this way, we arc led to the following triplet 
states: 
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The transition matrix elements induced by stag- 
gered spin operators H^\Pi) are expressed by the 
reduced matrix elements M^ n \l,P,0) (US) and 
M^ n '(2, P, 1) and appropriate Clebsch-Gordan co- 
efficients; the latter are absorbed in the definitions 
(|C10jl - (|Cl?P of the triplet states. Ap art f ro m ro- 
tations, there are two different states (|C1 1|) . (|C1G[) 
with three triplets. They turn out to be orthogonal 
to each other: 
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b) In this step we rotate the triplet states (|C1|) - ()C7I) 
and l(ClO)) - ((CT7)) as described in Ijg^}. In this way 
we obtain the decompositions 



d) In this step we compute the scalar products (I8.12p . 
The rt-independent 7x7 matrix (|8.3[) is given as 
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